CONTROL DEVICE FOR ELECTRIC MOTOR 



5 BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates to a control device for an electric motor. 
More in detail, this invention relates to a sensor-less control technique for a 
three-phase synchronous motor equipped with permanent magnets for a 
10 rotor, which will be called as "internal permanent magnet (IPM) motor" 
hereinafter. 

2. Description of Related Art 

In order to control the rotating speed of an electric motor to be driven by 
an inverter, it is necessary to detect the positions (phase) of poles of a rotor. 

15 As a method of detecting the position of a rotor without any sensor, there is 
an article titled "Evaluation for estimated accuracy in the positional 
sensor-less field detecting method for PM (permanent magnet) motor, using 
current vector locus 57 (No. 683, at National Convention in 1994 of Society 
for Industrial Application, Institute of Electric Engineers). The article 

20 describes a technique to detect a long axis of current vector locus in a 

high-frequency space by superimposing a high-frequency current on a drive 
current, thereby detecting the phase of a rotor. 

SUMMARY OF THE INVENTION 
25 The above-mentioned technique has a problem as follows. First of all, 

it is generally noted an electric motor has two characteristic features. 
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First, the more a q-axis current is increased, the more a q-axis 
inductance is reduced due to magnetic saturation, thereby causing a 
high-frequency current vector locus to approach a perfect circle. 

-- CD 

5 Second, due to the magnetic saturation, a position allowing magnetic 

flux to flow easiest and another position allowing magnetic flux to flow 
hardest move together with the phase of current. — (D 

With the delay in either errors or indication of a drive current as a start, 
if a q-axis current is increased due to the above characteristic feature (2), 

10 then a phase of detection is delayed, thereby causing a further increasing of 
the q-axis current. As a result, when an elliptical current becomes a perfect 
circle due to the above characteristic feature ® and when the correlation of 
the characteristic feature (2) approaches approximately 1 so that the long 
axis of an ellipse does not represent a d-axis but a current phase, the control 

15 accuracy is deteriorated at the motor's driving under a high load condition. 

The present invention is provided to improve the art described above. 
An object of the present invention is to provide a control device for an 
electric motor that enables a sensor-less operation even under a high load 
driving condition. 

20 According to the present invention, the above-mentioned object is 

accomplished by a control device for an electric motor, comprising: an 
inverter circuit connected to a direct-current power source to convert a 
power from the direct-current power source to an alternating current power, 
thereby driving the electric motor by the alternating current power; a control 

25 unit for controlling the operation of the inverter circuit on the basis of a 
target torque; a superimposing unit for superimposing a superimposed 
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current on a drive current for driving the electric motor, the superimposed 
current having a frequency different from a frequency of the drive current; a 
feature target generating unit for calculating a target value of a feature based 
on at least one of the length of a long axis of a current vector locus and the 
5 length of a short axis thereof, corresponding to the target torque; a 

separating unit for separating the superimposed current from a current 
supplied to the electric motor; an actual feature detecting unit detecting an 
actual value of the feature based on at least one of the length of a long axis 
of a current vector locus of the superimposed current separated by the 

10 separating unit and the length of a short axis thereof, the superimposed 
current having a different frequency from that of the drive current; and a 
phase angle detecting unit for detecting a phase angle of the electric motor 
on the basis of the target value calculated by the feature target generating 
unit and the actual value detected by the actual feature detecting unit, 

15 wherein the control unit controls the operation of the inverter circuit on the 
basis of the phase angle detected by the phase angle detecting unit. 

In the present invention, there is also provided a control device for an 
electric motor, comprising: an inverter circuit connected to a direct-current 
power source to convert a power from the direct-current power source to an 

20 alternating current power, thereby driving the electric motor by the 

alternating current power; a control unit for controlling the operation of the 
inverter circuit on the basis of a target torque; a feature target generating unit 
for calculating a target value of a feature based on at least one of the length 
of a long axis of a current vector locus and the length of a short axis thereof, 

25 corresponding to the target torque; a superimposing unit for superimposing a 
superimposed current on a drive current for driving the electric motor, the 
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superimposed current having a different frequency from that of the drive 
current; a separating unit for separating the superimposed current from a 
current supplied to the electric motor; a d-axis phase angle detecting unit for 
detecting a d-axis phase angle from a value of the superimposed current 
5 separated by the separating unit; an actual feature detecting unit for 

detecting an actual value of the feature based on at least one of the length of 
a long axis of a current vector locus of the superimposed current and the 
length of a short axis thereof; a correction angle calculating unit for 
calculating a correction angle on the basis of the target value calculated by 

10 the feature target generating unit and the actual value detected by the actual 
feature detecting unit; and a phase angle correcting unit for correcting the 
d-axis phase angle by the correction angle; wherein the control unit controls 
the operation of the inverter circuit on the basis of the d-axis phase angle 
corrected by the phase-angle correcting unit. 

15 Further, there is provided a method for controlling an electric motor 

having an inverter circuit connected to a direct-current power source to 
convert a power from the direct-current power source to an alternating 
current power for driving the electric motor, the method comprising: 
calculating a target value of a feature based on at least one of the length of a 

20 long axis of a current vector locus and the length of a short axis thereof, 
corresponding to a target torque for the electric motor; superimposing a 
superimposed current on a drive current for driving the electric motor 
thereby detecting an actual value of the feature based on at least one of the 
length of a long axis of a current vector locus of the superimposed current 

25 and the length of a short axis thereof, the superimposed current having a 

different frequency from that of the drive current; detecting a phase angle of 
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the electric motor on the basis of the target value and the actual value; and 
controlling the operation of the inverter circuit on the basis of the phase 
angle. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 is a block diagram showing the whole structure of a motor control 

system equipped with a control device in accordance with one embodiment 

of the present invention; 

Fig. 2 is a block diagram showing the details of the control device of Fig. 

1, in accordance with the first embodiment of the present invention; 

Fig. 3A is a diagram showing a voltage vector locus drawing a perfect 

circle, and Fig. 3B is a diagram showing a current vector locus drawing an 

ellipse; 

Figs. 4 A, 4B are flow charts of the process to start the torque control in 
case of no magnetic saturation (except a high load zone); 

Fig. 5 is a diagram to explain a method of detecting phases (£> i a and 
<M/3; 

Fig. 6 is a diagram showing a difference 6 e between a long axis of an 
ellipse and a d-axis at magnetic saturation in a current vector locus; 

Fig. 7 is a diagram showing a difference 6 e between a long axis of an 
ellipse and a d-axis at magnetic saturation in a voltage vector locus; 

Fig. 8 is a diagram showing the relationship between the length (a) of a 
long axis of an elliptic current as a characteristic feature and the current 
phase j3 ; 

Fig. 9 is a diagram showing the relationship between the length (a) of a 
long axis of an elliptic current as a characteristic feature and the current 



phase /3 ; 

Fig. 10 is a diagram showing the relationship between the length (a) of a 
long axis + the length (b) of a short axis of an elliptic current as a 
characteristic feature and the current phase /3 ; 

Fig. 11 is a diagram showing the relationship between the length (a) of a 
long axis + the length (b) of a short axis of an elliptic voltage as a 
characteristic feature and the current phase /3 ; 

Fig. 12 is a diagram showing the relationship between the length (b) of a 
short axis of an elliptic current as a characteristic feature and the current 
phase /3 ; 

Fig. 13 is a diagram showing the relationship between the length (a) of a 
short axis of an elliptic voltage as a characteristic feature and the current 
phase j3 ; 

Fig. 14 is a diagram showing the relationship between the length (b) of a 
short axis / the length (a) of a long axis of an elliptic current as a 
characteristic feature and the current phase j3 ; 

Fig. 15 is a diagram showing the relationship between the length (b) of a 
short axis / the length (a) of a long axis of an elliptic voltage as a 
characteristic feature and the current phase j3 ; 

Fig. 16 is a diagram showing the relationship between the length (a) of a 
long axis X the length (b) of a short axis of an elliptic current as a 
characteristic feature and the current phase j3 ; 

Fig. 17 is a diagram showing the relationship between the length (a) of a 
long axis X the length (b) of a short axis of an elliptic voltage as a 
characteristic feature and the current phase j3 ; 

Fig. 18 is a diagram regarding the length (a) of a long axis and the length 
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(b) of a short axis of an elliptic current as a characteristic feature, showing 
the relationship between 7~ (a 2 + b 2 ) and the current phase j3 ; 

Fig. 19 is a diagram regarding the length (a) of a long axis and the length 
(b) of a short axis of an elliptic voltage as a characteristic feature, showing 
5 the relationship between <J~ (a 2 + b 2 ) and the current phase /3 ; 

Fig. 20 is a diagram regarding the length (a) of a long axis and the length 
(b) of a short axis of an elliptic current as a characteristic feature, showing 
the relationship between <f~ (a 2 + b 2 ) (a + b) and the current phase j3 ; 

Fig. 21 is a diagram regarding the length (a) of a long axis and the length 
10 (b) of a short axis of an elliptic voltage as a characteristic feature, showing 
the relationship between J~ (a 2 + b 2 ) (a + b) and the current phase j3 ; 

Fig. 22 is a block diagram showing the details of the control device of 
Fig. 1, in accordance with the second embodiment of the present invention; 

Fig. 23 is a block diagram showing the details of the control device of 
15 Fig. 1, in accordance with the third embodiment of the present invention; 

Fig. 24 is a flow chart of the process to start the torque control in case of 
no magnetic saturation (except a high load zone) in the third embodiment; 

Fig. 25 is a diagram to explain a method of detecting phases <k> i OL and 
<t> i j3 in the third embodiment; 
20 Fig. 26 is a block diagram showing the details of the control device of 

Fig. 1, in accordance with the fourth embodiment of the present invention; 

Fig. 27 is a block diagram showing the details of the control device of 
Fig. 1, in accordance with the fifth embodiment of the present invention; 

Fig. 28 is a block diagram showing the details of the control device of 
25 Fig. 1, in accordance with the sixth embodiment of the present invention; 

Fig. 29 is a block diagram showing the details of the control device of 
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Fig. 1, in accordance with the seventh embodiment of the present invention; 

Fig. 30A is a diagram showing a current vector locus drawing a perfect 
circle, and Fig. 30B is a diagram showing a voltage vector locus drawing an 
ellipse; 

Fig. 31 is a flow chart of the process to start the torque control in case of 
no magnetic saturation (except a high load zone) in the seventh 
embodiment; 

Fig. 32 is a diagram to explain a method of detecting phases (f> v a and 
<t> v j3 in the seventh embodiment; 

Fig. 33 is a block diagram showing the details of the control device of 
Fig. 1, in accordance with the eighth embodiment of the present invention; 

Fig. 34 is a block diagram showing the details of the control device of 
Fig. 1, in accordance with the ninth embodiment of the present invention; 

Fig. 35 is a flow chart of the process to start the torque control in case of 
no magnetic saturation (except a high load zone) in the ninth embodiment; 

Fig. 36 is a diagram to explain a method of detecting phases <t> vd and 
<t> vq in the ninth embodiment; 

Fig. 37 is block diagram showing the details of the control device of Fig. 
1, in accordance with the tenth embodiment of the present invention; 

Fig. 38 is a block diagram showing the details of the control device of 
Fig. 1, in accordance with the eleventh embodiment of the present invention; 
and 

Fig. 39 is a block diagram showing the details of the control device of 
Fig. 1, in accordance with the twelfth embodiment of the present invention. 

DESCRIPTION OF THE PREFERRED EMBODIMENT 
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Referring to the accompanying drawings, embodiments of the present 
invention will be described below. 

Fig. 1 is a block diagram showing the whole structure of a motor control 
system, in accordance with one embodiment of the present invention. 

In Fig. 1, a control device 1 whose details will be described later is 
constructed so as to input respective signals of a current sensor 4 and a 
voltage sensor 6, further calculate a PWM (Pulse Width Modulation) 
command 7 and transmit it to an inverter circuit 2. The inverter circuit 2 
converts a direct-current power of a power source 5 to a three-phase power 
corresponding to the PWM command 7 thereby to drive an IPM (Internal 
Permanent Magnet) motor 3. The current sensor 4 detects two-phase 
currents (e.g. U-phase current and V-phase current) of the three-phase power 
transmitted from the inverter circuit 2 to the IPM motor 3. While, the 
voltage sensor 6 detects an output voltage of the power source 5, that is, an 
input voltage inputted to the inverter 2. Detection values of the current 
sensor 4 and the voltage sensor 6 are transmitted to the control device 1 to 
calculate the PWM command 7. Noted that with the establishment of a 
relationship (U+V+W=0) about the three-phase currents, if only two 
currents are detected, the remaining current can be calculated automatically. 

Fig. 2 is a block diagram showing the details of the control device 1, in 
accordance with the first embodiment of the invention. 

In Fig. 2, a current/characteristic feature target generating unit 10 
generates a target value in d-axis current, a target value in q-axis current and 
a characterizing amount based on a target torque (e.g. manipulating amount 
for accelerator pedal) provided from outside. Noted that "characteristic 
feature" is a value based on at least one of the length of a long axis and the 



length of a short axis of a current vector locus. As the above value, for 
example, there are the length (a) of a long axis, a length (b) of a short axis, a 
+ b, aXb, ^( a 2 + b 2 ), a/b, V^(a 2 + b 2 )^(a + b), etc. which will be 
described later in detail. 

A d-axis/q-axis converting unit 11 calculates a d-axis current and a 
q-axis current by the three-phase currents obtained by the current sensor 4 of 
Fig. 1 and a phase angle 0 (its details described later). 

A current control unit 12 (conversion of current to voltage) carries out, 
for example, a PI control (known proportional/integral control) to calculate a 
d-axis voltage command to allow the target value in q-axis current to 
coincide with a d-axis current and a q-axis voltage command to allow the 
target value in q-axis current to coincide with a q-axis current. 

Note, the above current control may contain a noninterference control 
and a low-pass filter to remove a high-frequency component in the 
three-phase currents. 

A three-phase converting unit 13 calculates three-phase voltages 
commands on the basis of the d-axis voltage command, the q-axis voltage 
command and a phase angle 6 and further outputs the so-obtained 
three-phase voltages commands. 

A high-frequency rotating voltage-generating unit 14 (superimposing 
means for superimposing a superimposed current of higher frequency than 
that of a drive current) generates a rotating voltage of high frequency to 
transmit it to a PWM command-generating unit 15. Here noted that the 
above "rotating voltage" means a three-phase equilibrium voltage, 
representing a state that its voltage vector locus becomes a perfect circle 
when converted to an a / j3 -coordinate system as two-phase AC 
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(alternating current) coordinate system, as shown in Fig. 3A. Since this 
rotating voltage is a high frequency voltage that does not synchronize the 
IPM motor 3, there is no possibility that the IPM motor 3 is rotated by the 
rotating voltage. 

While, if there is a difference in inductance between flux current and 
torque current as the IPM motor, its current vector locus becomes an ellipse 
expanding to a direction of d-axis, as shown in Fig. 3B. 

The PWM command generating unit 15 inputs a voltage obtained by 
superimposing the rotating voltage from the high-frequency rotating 
voltage-generating unit 14 on the three-phase voltage command and a DC 
voltage brought by the voltage sensor 6 of Fig. 1, thereby to generate the 
PWM command 7. The inverter circuit 2 of Fig. 1 is controlled by the 
PWM command 7, while the IPM motor 3 is driven by a three-phase power 
generated from the inverter circuit 2. 

A frequency-separating unit 16 (means for separating a superimposed 
current) separates a high-frequency current from a three-phase current of the 
current sensor 4 and further transmits the high-frequency current to a 
3/2-phase converting unit 17. 

The 3/2-phase converting unit 17 converts a three-phase current of high 
frequency from the frequency separating unit 16 to a two-phase current in 
the a 1 j3 - coordinate system of Fig. 3B and further transmits the 
so-converted two-phase current in the a/ (3 - coordinate system to a d-axis 
detecting unit 18. 

Regarding a -axis component i a and j3 -axis component i j3 of the 
current vector as a result of converting to the a / j3 - coordinate system by 
the 3/2-phase converting unit 17, the d-axis detecting unit 18 (d-axis phase 



angle detecting unit) detects respective amplitudes I CL , I $ by detecting their 
peaks of the components i a , i 0 and further detects respective phases <$> i 
Q! , 0 i /3 of the superimposed currents by detecting their zero-cross 
periods of the components i a , i j3 , as shown in Fig. 5. In the a / /3 - 
coordinate system of Fig. 3B, an angle 6 from the a -axis to the long axis 
(d-axis) of an ellipse is obtained by the following expressions (1). 



0= - tan ] (-^-), (A<B) 
2 y A-B } - v ' 



6=1 tan 1 (-|^-)+(7t/2),(A>B,H>0) 



<9 = 1 tan 1 (-^-)-(7C/2),(A>B,H<0) 



A-B 



6 = 
6 = 



71/4, 
71/4, 



(A=B, H>0) 
(A=B, H<0) 



Note: A=l/( la) 2 , B=l/( I j8 f, H= 

ia=IO!sin(cot+^ ia ) 
i j3 = 1/3 sin (a>t+T/> i/3 ) 



- (1) 

__ cosQ/> ia-il)i(3) 



IaxIjS 



By the above expression (1), an angle 0 between the a -axis and the 
/3 -axis is obtained within a zone of ±90 degrees. 

In order to expand this angle 6 up to a zone of ± 180 degrees, there is 
employed a premise that there is no change in phase exceeding 90 degrees in 
accordance with a continuous detection. That is ? as the later-mentioned 
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flow chart of Fig. 4 starts, for example, every 100 U m, a phase does not 
change more than 90 degrees generally. Accordingly, in the continuous 
detection, when a phase changes from 89 degrees to — 89 degrees, it is 
assumed that the phase changes from 89 degrees to 91 degrees (= — 89 + 180 
degrees). Similarly, when a phase changes from — 89 degrees to 89 
degrees, it is assumed that the phase changes from — 89 degrees to — 91 
degrees (=89 + 180 — 360 degrees). Consequently, it is possible to detect 
an angle 6 between a -axis and d-axis within a zone of ± 180 degrees. 

Noted that the initial phase at the first detection cannot be determined 
because of no presence of comparative object in change. In detail, when 
the calculation at the first detection results in 30 degrees, it is unclear 
whether the initial phase is 30 degrees or — 150 degrees. Therefore, it is 
established at the first detection that the angle 6 zones from — 90 degrees 
to 90 degrees temporarily, in the above case, 30 degrees. Subsequently, at 
the later-mentioned polar judgment of N/S, the angle 0 is detected within 
the zone of +180 degrees. 

An angle 0 from the a -axis to the long axis (d-axis) of an ellipse 
detected at the d-axis detecting unit 18 is established as a d-axis phase angle 
0 0 . 

A characteristic-feature calculating unit 19 ("characteristic feature" 
actual value detecting means) calculates an actual characteristic feature (i.e. 
parameter) by both of the a -axis current and the j3 -axis current, of which 
details will be described later. 

A characteristic-feature control unit 20 (correction angle calculating 
means) calculates a correction angle 6 ' by the actual characteristic feature 
and the above target value in the characteristic feature. That is, when the 
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actual characteristic feature is larger than the target value in the 
characterizing amount, it is carried out to advance the detection phase. 
While, when the actual characteristic feature is smaller than the target value, 
it is carried out to delay the detection phase. In this feedback control, PI 
control method is used to establish a control gain by way of experiment. 

For example, when the phase angle /3 of the drive current increases 
and the value of the characteristic feature decreases, the calculation of the 
correction angle 8 9 by the PI control method is carried out by using the 
following expression (2). 

Noted, similarly to the case of the characteristic feature of S~ (a 2 + b 2 ) 
(a + b), when the value of the characteristic feature increases together 
with the increasing of the phase angle j3 of the drive current, the 
expression (2) is changed to the expression (2') described later. 



0'(s)= -(K P + ^){ft*(s) - ft(s)} - (2) 



8 ' : correction angle 

ft* : target value in characteristic feature 
ft : characteristic feature 
K p : proportional gain 
Kj : integral gain 
s : Laplace operator 



6"(s) = (K p+ ^){ft*(s) - ft(s)} - (2') 

3 
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As for the target value in characterizing amount at the 
current/characteristic feature target generating unit 10, it is noted that there 
are methods of changing the target value together with the target value in 
current, making the target value constant and so on. In any case, since the 
correction by the correction angle 9 ' is not applicable to a low load zone 
of small current, the activation/inactivation of control is executed depending 
on the amplitude in current. 

That is, the switch SW1 is turned ON only when the IPM motor 3 is 
driven under a high load zone where its drive current detected by the current 
sensor 4 is more than a predetermined value. 

When the switch SW1 is turned ON, the d-axis phase angle 9 0 from 
the d-axis detecting unit 18 is added to the correction angle 9 ' from the 
characteristic-feature control unit 20 at an accumulator 26 (phase-angle 
correcting unit), whereby the d-axis phase angle 9 0 is corrected with the 
correction angle 9 providing a phase angle 9 . Then, the phase angle 
9 after correction is transmitted to the d-axis/q-axis converting unit 11 and 
the three-phase converting unit 13 for the above-mentioned calculation. 

In the circuit of Fig. 2, therefore, when the IPM motor 3 is driven under 
the high load zone where its drive current detected by the current sensor 4 is 
more than the predetermined value, the phase angle 9 corrected on the 
basis of the characterizing amount is used for calculation. When the IPM 
motor 3 operates except the high load zone, the d-axis phase angle 9 0 
generated from the d-axis detecting unit 18 is used as it is. 

As mentioned above, if only constructed so as to correct the d-axis phase 
angle 9 0 detected by the d-axis detecting unit 18 only at a high load zone 
where the drive current for the motor is more than a predetermined value, it 
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is possible to stop the calculations of the correction angle 6 9 at the 
characteristic-feature calculating unit 19 and the characteristic-feature 
control unit 20 in such a low load zone that requires no correction for the 
phase angle, allowing the calculation load for the control device to be 
5 reduced. 

Fig. 4 is a flow chart of the process to start the torque control in case of 

no magnetic saturation (out of a high load zone). 

First, in Fig. 4A, at step SI, it is executed to detect a temporary d-axis 

by passing a high-frequency rotating current corresponding to a 
10 high-frequency rotating voltage from the high-frequency rotating voltage 

generating unit 14. 

At step S2, it is executed to detect a true d-axis by the N/S polar 

judgment for specifying N-pole and S-pole of magnetic poles. As for the 

N/S polar judgment, there are a variety of methods of: (1) applying such 
15 high-frequency voltage as a large voltage to cause a saturation of flux 

thereby to detect a direction to deviate from the center of an ellipse; (2) 

applying a low-frequency sinusoidal voltage to only a d-axis at stopping 

thereby to compare a positive magnitude of the d-axis current with its 

negative magnitude; (3) applying positive and negative currents to only a 
20 d-axis thereby to compare amplitudes of high-frequency currents in both 

cases with each other, and so on. In this embodiment, any one of these 

methods is employed. 

On the above detection of a true d-axis, namely, d-axis phase angle, the 

torque control is performed with use of the true d-axis at step S3. The 
25 contents of the torque control have been described on and after the previous 

descriptions about the current control unit 12. 



16 



Next, at step S4 of Fig. 4B, it is executed to separate a high-frequency 
component from the current detected by the current sensor 4 (by the 
frequency separating unit 16). Then, at step S5, the high-frequency current 
is converted to the a I j3 -coordinate system to realize the state of Fig. 3B 
(by the 3/2-phase converting unit 17). 

At step S6, as for both a -axis component i a and /3 -axis component 
i j3 of the current vector converted to the a/ j3 - coordinate system, it is 
executed to detect respective amplitudes I OL , I j3 by detecting their peaks of 
the components i a , i j3 and also executed to detect the superimposed 
current phases <f> i a , <t> i j3 by detecting their zero-cross periods of the 
components i 05 , i j3 , as shown in Fig. 5. 

At step S7, it is executed to calculate the initial angle 0 within the 
zone of ±90 degrees provisionally. Noted that the calculation result at 
this time is represented by an angle 6 2 while the calculation result at the 
previous routine is represented by an angle 6 lm 

At steps S8 to S10, it is executed to expand the initial angle 0 within 
the zone of ± 90 degrees to an angle within the zone of ± 180 degrees. In 
order to expand this angle 9 up to the zone of +180 degrees, there is 
employed a premise that there is no change in phase exceeding 90 degrees in 
accordance with a continuous detection. That is, in accordance with the 
continuous detection, when a phase changes from 89 degrees to — 89 
degrees, it is assumed that the phase changes from 89 degrees to 91 degrees 
(= — 89 + 180 degrees). Similarly, when a phase changes from — 89 
degrees to 89 degrees, it is assumed that the phase changes from — 89 
degrees to — 91 degrees (=89 + 180 — 360 degrees). Consequently, it is 
possible to detect an angle 6 within a zone of ±180 degrees. However 
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noted that the initial phase at the first detection cannot be determined. In 
detail, if the calculation at the first detection results in 30 degrees, it is 
unclear whether the initial phase is 30 degrees or — 150 degrees. 
Therefore, it is established at the first detection that the angle 6 zones from 
— 90 degrees to 90 degrees temporarily. 

At step S8, it is judged whether the present calculation is the initial 
calculation or not. If the judgment at step S8 is Yes (the initial calculation), 
then the routine goes to step Sll to output the calculation result (value) as it 
is. On the other hand, if the judgment at step S8 is No (not the initial 
calculation), then the routine goes to step S9. 

At step S9, it is judged whether the absolute value of a difference 
between the previous calculation value (angle) 6 i and the present 
calculation value (angle) 6 2 is more than 90 degrees. If the judgment at 
step S9 is No (less than 90 degrees), then the routine goes to step Sll to 
output the present calculation value (angle) 6 2 as a detected phase angle 6 0 . 
If the judgment at step S9 is Yes (more than 90 degrees), then the routine 
goes to step Sll to output a value obtained by adding 180 degrees to the 
present calculation value (angle) 9 2 as a detected phase angle 6 0 . 

The above-mentioned calculation (from step S6 to step Sll) is 
performed by the d-axis detecting unit 18. 

We now describe the correction of the phase angle by means of the 
characterizing amount as the feature of the present invention. 

The correction of the phase angle at the d-axis detecting unit 18 and the 
characteristic-feature control unit 20 is carried out when the magnetic 
saturation arises (at high-load zone). 

When the magnetic saturation arises, a difference 0 e between the long 
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axis of an ellipse and the d-axis becomes as shown in Fig. 6. Here noted 
that a term "ellipse" means an ellipse (elliptic current) of the current vector 
locus shown in Fig. 3B. In Fig. 6, a longitudinal axis designates the 
difference 6 e between the long axis and the d-axis, a transverse axis /3 
designates a phase angle (i.e. phase angle from the q-axis) of the drive 
current and sequential lines denote the characteristics of various drive 
currents i a =S~ (i d 2 + i q 2 ). When the drive current value is small, there is no 
change in a fact that the long axis of the ellipse substantially represents the 
d-axis in spite of a little error therebetween (the zone from about 100A to 
300A in Fig. 6). While, when the drive current value is large, the long axis 
of the ellipse moves in accordance with the current phase angle, enhancing 
the correlation with the current phase angle apart from the correlation with 
the d-axis. In this situation, it is impossible to detect a current phase. 
Therefore, in case of magnetic saturation (at high load zone), the other 
characteristic feature obtained by an elliptic current is employed. 

Fig. 7 shows the characteristics similar to that of Fig. 6, with respect to a 
voltage vector locus (shown in Fig. 30B). As shown in Fig. 30B, since the 
position of (the long axis + 90 degrees) corresponds to the d-axis in the 
voltage vector locus (elliptic voltage), Fig. 7 exhibits the difference 6 e 
between (the long axis + 90 degrees) and the d-axis. 

Figs. 8, 10, 12, 14, 16, 18 and 20 are views showing the relationships 
between various characteristic features and the drive current phase angle 13 
under the similar condition as that of Fig. 6. In these figures, Fig. 8 shows 
the characteristic feature of the length a of the long axis in an elliptic current, 
Fig. 10 the characteristic feature of the total length of (a + b), Fig. 12 the 
characteristic feature of the length b of the short axis, Fig. 14 the 
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characteristic feature of the calculation of (b/a), Fig. 16 the characteristic 
feature of the calculation of (a X b), Fig. 18 the characteristic feature of 
the calculation of ^T(a 2 + b 2 ), and Fig. 20 shows the characteristic feature of 
the calculation of ,T(a 2 + b 2 )-^(a + b). Similarly, Figs. 9, 11, 13, 15, 17, 
19 and 21 are views showing the relationships between various 
characteristic features and the drive current phase angle /3 under the 
similar condition as that of Fig. 7. In these figures, Fig. 9 shows the 
characteristic feature of the length a of the long axis in an elliptic voltage, 
Fig. 11 the characteristic feature of the total length of (a + b), Fig. 13 the 
characteristic feature of the length b of the short axis of an elliptic voltage, 
Fig. 15 the characteristic feature of the calculation (b/a), Fig. 17 the 
characteristic feature of the calculation (a X b), Fig. 19 the characteristic 
feature of the calculation of 4~ (a 2 + b 2 ), and Fig. 21 shows the characteristic 
feature of the calculation of <T (a 2 + b 2 ) -5- (a + b). Note, embodiments 
adopting the above voltage vector locus will be described in Fig. 29, 
thereinafter. 

The above characteristic-feature control unit 20 is provided to control 
these characteristic features to respective predetermined values, i.e. 
respective targets determined at the current/characteristic feature target 
generating unit 10. In the characteristic-feature control unit 20, a 
correction angle 9 ' for phase is calculated to correct a d-axis phase angle 
6 o by means of PI control etc. 

In accordance with the properties of Figs. 8, 10, 12, 14, 16, 18 and 20, 
the characteristic features have downward-sloping features as the drive 
current phase angle j3 increases. Therefore, if an actual characteristic 
feature is larger than its target value, the detection phase is advanced 
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(increasing of the correction angle 9 ') so as to advance the rotating 
coordinate system while fixing the drive current phase angle j3 at a 
constant, thereby producing an effect to put the current phase angle j3 
forwardly. While, if an actual characteristic feature is smaller than its 
5 target value, the detection phase is delayed (decreasing of the correction 
angle 9 ') so as to retard the rotating coordinate system while fixing the 
drive current phase angle j3 at a constant, thereby producing an effect to 
put the current phase angle /3 backwardly. Then, the control gains are 
determined by way of experiments. Note, as the feature of Fig. 20 

10 (regarding the characteristic feature of the calculation of 7" (a 2 + b 2 )-^ (a + 
b)) ? if a characteristic feature has a feature of upward-sloping as the drive 
current phase angle # increases, the above-mentioned relationship 
between advance and delay is reversed. 

However, the correction by the correction angle 9 ' is carried out only 

15 in such a high load zone as to cause the magnetic saturation. In other 
words, the switch SW1 is turned ON when the drive current for the IPM 
motor 3 is more than a predetermined value, and also turned OFF when the 
drive current is less than the predetermined value. 

Although there are supposed length of long axis (a), length of short axis 

20 (b), a + b, </"" (a 2 + b 2 ), etc. as the characteristic features obtained by the 

magnitude of an elliptic current, the using of characteristic features (e.g. b/a ? 
*T (a 2 + b 2 )-r-(a + b), etc.) with no unit (i.e. ratio) enables more precise 
control. 

In adopting the characteristic feature of the calculation (b/a), it is 
25 possible to prevent an elliptic current from being of a perfect circle. While, 
in case of adopting the other characteristic features, their target values have 



21 



to be established in order to prevent a circular elliptic current. 

In case of the characteristic feature of (aXb), as shown in Fig. 16, since 
the current phase angle at a maximum torque zones from 30 to 50 degrees 
approximately, the inclinations for every currents are so linear as to facilitate 
the control of the characteristic feature in spite of the presence of minor 
differences in inclination. This characteristic feature is applicable to the 
characteristic features of a + b and <T (a 2 + b 2 ), as well. However, as the 
characteristic feature depends on the motor, not all motor exhibits preferable 
characteristic feature. Noted, in case of the characteristic feature of (a X b), 
it has an advantage of reducing the calculation load because of no 
calculation of the square root (V~ ). 

The above-mentioned characteristic features are obtained by the 
following expressions. Thus, providing that alphabet (a) represents the 
length of a long axis of an ellipse and alphabet (b) represents the length of a 
short axis of the ellipse, both of the lengths (a) and (b) are calculated by the 
following expression (3). 

a= \\T {I, 2 +I* 2 +2I 0 I fl I sin(^ ia -^i, ) I } 

+ AT{Ia 2 +I 0 2 - 2I 0 I fl I sin(i/> ia -4> i/3 )|}] 
b= \[f {I a 2 +I e 2 +2I a l/ I sinO ia -i/> i/3 )l} 

-V"{Ia 2 +I/3 2 - 2I a I 3 | Sin(^ ia -7/> i/3 ) I }] 

- (3) 



As the target value, the current/characteristic feature target generating 
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unit 10 includes the characteristic features obtained by way of experiments 
in advance, in the form of a table for every torque. Similarly to the target 
values in current, the same unit 10 searches and calculates the characteristic 
feature corresponding to the target value in torque. Besides, since the 
actual current is different from the target value in current, it is desirable to 
apply a low-pass filter considering the responsibility in current control to the 
target value for the characteristic feature. Alternatively, the unit 10 may be 
adapted so as to calculate a target value for the characteristic feature from an 
actual current while using a table for every current instead of the above table 
for every torque. 

Now, we provide respective expressions for the above-mentioned 
calculations, together. 
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(4) 



id , i q : discrete values of d-axis and q-axis currents [A] 
111 3 lv ? lw : discrete values of three-phase currents [A] 

(note: the member "i w " may be omissible on the relationship of '% 
iu-iv ") 

0 : phase generated from the d-axis detection unit 

v ; (s) = (Kp d + K id /s){ - (s) - id (s)} 

y ; (s) = (K pq + K iq /s){ j; (s) - iq (s)} 

- (5) 
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Id 



Kjd, Kj q 



: command values of d-axis and q-axis currents [A] 
: command values of d-axis and q-axis voltages [V] 
: Laplace operator 

: proportional gains of d-axis and q-axis 
: integral gains of d-axis and q-axis 
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(6) 



Vd ' V, 



Vd' v q 
outputs [V] 

15 i d , i q 

CO 



command values of d-axis and q-axis voltages [V] 
command values of d-axis and q-axis voltages as PI control 

values of d-axis and q-axis currents [A] 
inductance of d-axis and q-axis currents [H] 
induced voltage constant [Wb] 
angular speed (electrical degree) [rad/s] 
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•(7) 



y* , y* , y* i command values of three-phase voltages [V] 
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- (8) 

v dc : discrete value of DC voltage [V] 

t' t ' t» : three-phase PWM duty widths (ON-period for upper arm) 
T 0 : PWM cycle (e.g. 100 [ fi s] at 10 [kHz]) 



a + b= f {I 0 2 +I e 2 +2I a I e | sin(7/> ia -i/> i(3 ) | } 

... (9) 

^(a 2 + b 2 )= 4~ (I a 2 +I, 2 ) 

••• (10) 

ab = I a I 0 | sin(T/>i tt -T/» i/5 ) I 

- (11) 



As mentioned above, according to the first embodiment of the invention, 
since the phase angle is corrected by the characteristic feature based on at 
least one of the length of a long axis and the length of a short axis of a 
current vector locus, it is possible to detect the phase angle without any 
sensor even when the motor is driven within a high load zone causing 
magnetic saturation. 

In addition, it is also possible to absorb a delay in detecting the phase 
due to a change in rotating speed of the motor. 

In case of employing the product (aXb) as the characteristic feature, it 
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has advantages of improving the characteristic feature in the vicinity (about 
30 to 60 degrees) of current phase angle for generating the maximum torque 
as similar to the other characteristic features (a + b, J~ (a 2 + b 2 )) and 
reducing the calculation load. 

Next, Fig. 22 is a block diagram, showing the details of the control 
device 1 of Fig. 1, in accordance with the second embodiment of the 
invention. 

In this embodiment, when the motor is driven in a high load zone 
causing the magnetic saturation, only the characteristic-feature control unit 
20' detects a phase of the motor, while a d-axis phase angle 6 0 from the 
d-axis detecting unit 18 is not used. 

In this case, the characteristic-feature control unit 20' outputs an angular 
speed co '[rad/s]. Then, the angular speed co ' is integrated to calculate a 
phase angle 6 ' [rad] by an integrating circuit 21. That is, the 
characteristic-feature control unit 20' calculates an angular speed co '[rad/s] 
and corrects a phase angle by means of PI control method etc. Since the 
characteristic quantity has a feature of downward inclination to the right due 
to the current phase £ , if an actual characteristic feature is larger than its 
target value, the detection phase is advanced (increasing of the angular 
speed co ') so as to advance the rotating coordinate system while fixing the 
at a constant, thereby producing an effect to advance the drive current phase 
angle /3 . While, if an actual characteristic feature is smaller than its target 
value, the detection phase is delayed (decreasing of the angular speed co ') 
so as to retard the rotating coordinate system while fixing the at a constant, 
thereby producing an effect to delay the drive current phase angle # . 
Then, the control gains are determined by way of experiments. 
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The calculation of the angular speed co ' is carried out by using the 
following expression (12), for example. On the other hand, as the feature 
of Fig. 20 (i.e. the characteristic feature of the calculation J~ (a 2 + b 2 )-^(a + 
b)), if a characteristic feature has a feature of upward-sloping as the drive 
5 current phase angle j3 increases, the angular speed co ' is obtained by the 
expression (12'). 

co(s) = -(K p+ ^){ft*(s) - ft(s)}- (12) 

co : angular speed (electrical angle) [rad/s] 
10 ft* : target value in characteristic feature 

ft : characteristic feature 
K p : proportional gain 
Kj : integral gain 
s : Laplace operator 

15 

co(s) = (K p + ^){ft*(s) - ft(s)} - (12') 

In this case, the relationship between the magnitudes of the actual 
characteristic feature/feature target value and the increase/decrease of the 
20 angular speed co is opposite to that of the case of having a feature of 
downward-sloping. 

With the exception of high load zone causing magnetic saturation, the 
switch SW1' is connected with the d-axis detecting unit 18 to transmit a 
d-axis phase angle 6 0 to the following circuit, in the form of a phase angle 
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9 . In the high load zone causing magnetic saturation, the switch SW1 is 
connected to the integrating circuit 21 to transmit a value as a result of 
integrating an angular speed co ' to the following circuit, in the form of a 
phase angle 9 \ 

The constitution of Fig. 22 has an advantage of eliminating attention that 
an elliptic current might be a perfect circle. 

Fig. 23 is a block diagram, showing the details of the control device 1 of 
Fig. 1, in accordance with the third embodiment of the invention. 

In this embodiment, an angular speed co is employed in the calculation 
of a d-axis phase angle 9 0 at the d-axis detecting unit 18. 

In Fig. 23, the frequency separating unit 16 inputs a d-axis current and a 
q-axis current after conversion by the d-axis/q-axis converting unit 11 and 
further separates high-frequency components (high-frequency d-axis/q-axis 
currents) from the above currents. 

At the d-axis detecting unit 18 both amplitude and phase of the current 
vector are determined by detecting respective peak values and zero-cross 
periods of the d-axis component id and the q-axis component iq of the 
current vector (after separating high-frequency components) converted onto 
the d/q-axes. Although an angle 6 between the d-axis and the long axis 
of an ellipse is similar to that of the above expression (1) expect to alter the 
members of I a , 1 /3 , <t>ia and <£i/3 to the members of Id, Iq, <$> id and 
<t> iq respectively, it is executed to multiply the angle 9 by a regular 
coefficient (i.e. proportional calculation) thereby outputting an angular 
speed co . Then, in the integrating circuit 22, the angular speed co is 
integrated to generate a d-axis phase angle 9 0 for output. 

As mentioned above, the angular speed co may be provided by the 
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proportional calculation while establishing the angle 9 from the d-axis to 
the long axis of an ellipse, of zero. Alternatively, the angular speed co 
may be provided by a PI (proportional integrating) calculation. Or again, 
providing that the unit 18 ' includes map-tables with every torque each 
defining the relationship between a torque detected by an experiment in 
advance and the corresponding angle 6 from the d-axis to the long axis of 
an ellipse, if the unit 18' is adapted so as to determine the angle 6 
corresponding to a target torque by means of the specified map-table, then it 
is possible to improve the accuracy of an estimated phase. 

The characteristic-feature control unit 20 outputs a correction angle 6 ' , 
as similar to Fig. 2. In a high-load zone where the switch SW1 is turned 
ON, the characteristic-feature control unit 20 corrects the d-axis phase angle 
6 o by the correction angle 8 ' to provide a phase angle 6 , as similar to 
Fig. 2. 

The detection of the long axis on the d/q-axes, which is embodied by the 
arrangement of Fig. 23, has an advantage of eliminating detection errors 
since the locus of a high-frequency current vector comes to a standstill while 
drawing an ellipse on the d/q-axes at regular revolutions of the motor. 

Fig. 24 is a flow chart of the process to start the torque control in case of 
no magnetic saturation (out of a high load zone), corresponding to the 
operation of Fig. 4B. Note, the operation shown in Fig. 24 follows the 
operation shown in Fig. 4A. 

In the flow chart of Fig. 24, at step S241, a current detected by the 
current sensor 4 is converted onto the d/q-axes by the d-axis/q-axis 
converting unit 11 and thereafter, the high-frequency component is separated 
from the current (frequency separating unit 16) thereby to establish a state as 
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shown in Fig. 3B. At sequent step S242, as for the d-axis component id 
and the q-axis component iq of the current vector converted in the d-q 
coordinate system, their amplitudes Id, Iq can be obtained by detecting their 
respective peak values of the components id, iq as shown in Fig. 25. 
Further, their superimposed current phases <t> id , <t> iq can be obtained by 
detecting their zero-cross times of the components id, iq. At sequent step 
S243, a phase angle 9 is calculated from d-axis to a long axis of an ellipse 
within ± 90 degrees. At step S244, P-control or Pi-control is performed in 
order to make the detected phase angle 9 coincident with the target phase 
angle 9 to output an angular speed co 0 . As the target phase angle 9 , a 
phase angle 0 degrees may be used, or the other value corresponding the 
torque that is experimentally obtained may be used in order to improve a 
phase angle estimation accuracy. At step S245, the angular speed oo 0 is 
integrated by the integration circuit to output the d-axis phase angle 9 0 . 
These operations (steps S243 to S244) are performed in the d-axis detecting 
unit 18'. 

In this embodiment to detect the long axis on the d/q-axes, however, the 
expressions (3), (9), (10) and (11) in the first and second embodiments will 
be modified to the following expressions (13), (14), (15) and (16), 
respectively. 

a= \[f {I d 2 +I q 2 +2I d I q | sin(7/> id --0 iq ) I } 

4- /~{I d 2 +I q 2 - 2I d I q I sin(^ id ~^ iq ) I }] 
b= {I d 2 +I q 2 +2I d I q 2 I sin^w-^iq) I } 
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-/"{I d 2 +I q 2 - 2I d I q | sin(T/> w — -0U,) I }] 

- (13) 



a + b= -f {I d 2 +I q 2 +2I d I q I sin(7/> id -^ iq ) | } 

■•• (14) 

f(a 2 + b 2 )= f (I d 2 +I q 2 ) 

- (15) 

ab = I d I q | sin(t/> id -i/> iq ) | 

(16) 

As the method of applying a high-frequency voltage, there is also 
another method of applying the voltage on the d/q-axes (in the rotating 
coordinate system) besides the application of the voltage through a 
three-phase converting unit (in the static coordinate system). That is, in the 
first to third embodiments, a high-frequency current is superimposed on a 
drive current by superimposing a high-frequency voltage command 
generated from the high-frequency rotating voltage generating unit 14 on a 
three-phase voltage command outputted from the three-phase converting 
unit 13. While, according to the above-mentioned method, it may be 
carried out to superimpose the high-frequency voltage command generated 
from the high-frequency rotating voltage generating unit 14 on a d/q-axes 
voltage command generated from the current control unit 12. Since the 
so-superimposed "high-frequency current" has a large frequency in 
comparison with the rotational frequency (i.e. frequency of the motor drive 
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current), the superimposition of a voltage command on the d/q-axes allows a 
slippage of frequency to be eliminated during the rotation of a motor 
although the superimposition of a high-frequency voltage command on a 
three-phase voltage command can exhibit little influence on the drive of a 
motor. 

Next, Fig. 26 is a block diagram, showing the details of the control 
device 1 of Fig. 1, in accordance with the fourth embodiment of the 
invention. In this embodiment, as similar to the above-mentioned third 
embodiment, the d-axis detecting unit 18' detects both amplitude and phase 
by detecting respective peak values and zero-cross periods of the d-axis 
component id and the q-axis component iq of the current vector (after 
separating high-frequency components) converted onto the d/q-axes. 
Further, at the d-axis detecting unit 18', it is executed to multiply the 
so-detected amplitude and phase by a regular thereby outputting an angular 
speed co . Then, in the integrating circuit 22, the angular speed co is 
integrated to generate a d-axis phase angle 0 0 for output. 

On the other hand, as similar to the second embodiment, the 
characteristic-feature control unit 20' outputs an angular speed co '[rad/s]. 
Then, the angular speed co ' is integrated to calculate a phase angle 6 ' 
[rad] by an integrating circuit 21. 

With the exception of high load zone causing magnetic saturation, the 
switch SW1' is connected with the integrating circuit 22 to transmit a d-axis 
phase angle 6 0 from the circuit 22 to the following circuit, in the form of 
a phase angle 6 . In the high load zone causing magnetic saturation, the 
switch SW1' is connected to the integrating circuit 21 to transmit a value 6 ', 
which is obtained as a result of integrating an angular speed co ', to the 



32 



following circuit, in the form of a phase angle 6 . 

Next, Fig. 27 is a block diagram, showing the details of the control 
device 1 of Fig. 1, in accordance with the fifth embodiment of the invention. 

In Fig. 27, the angular speed co from the d-axis detecting unit 18' is 
transmitted to the switch SW1' directly, while the characteristic-feature 
control unit 20' outputs an angular speed co ', as similar to the embodiment 
of Fig. 26. Then, the switch SW1 ' is connected with the 
characteristic-feature control unit 20' in a high-load zone causing magnetic 
saturation. Except the high-load zone, the switch SW1 ' is connected with 
the d-axis detecting unit 18'. Therefore, in the high-load zone, the 
integrating circuit 23 integrates the angular speed co ' from the 
characteristic-feature control unit 20' to output a phase angle 6 . Except 
the high-load zone, the integrating circuit 23 integrates the angular speed co 
from the d-axis detecting unit 18' to output a phase angle 6 . 

Note, in switching an angle 6 as shown in the embodiments of Figs. 22 
and 23, if there is a difference between the phase calculated at the d-axis 
detecting unit 18' and the phase calculated at the characteristic-feature 
control unit 20', there is produced a skip in phase. While, in case of 
switching an angular speed co in this embodiment, it is possible to solve 
the above problem of phase-skipping. 

Next, Fig. 28 is a block diagram, showing the details of the control 
device 1 of Fig. 1, in accordance with the sixth embodiment of the 
invention. 

In Fig. 28, the angular speed co outputted from the d-axis detecting 
unit 18' is transmitted to the accumulator 26 directly. While the 
characteristic-feature control unit 20" outputs a correction angular speed 
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co ". In a high-load zone where the switch SW1' is turned on, the 
correction angular speed co " is added to the angular speed co at the 
accumulator 26. 

The integrating circuit 23 integrates an output from the accumulator 26 
and further outputs a phase angle 9 . Therefore, except the high-load zone, 
the integrating circuit 23 integrates the angular speed co outputted from the 
d-axis detecting unit 18' to output the phase angle 9 . While, in the 
high-load zone, the integrating circuit 23 integrates a value obtained by 
correcting the angular speed co by the correction angular speed co " to 
output a phase angle 9 . Note, although the arithmetic expression of a 
correction angular speed co " is similar to the expressions (12) and (12'), as 
evidenced by these expressions (12) and (12'), the characteristic-feature 
control unit 20" operates in a manner that the target value in characteristic 
feature ft* accords with the characteristic feature ft, so that an output of the 
unit 20" becomes the correction angular speed co ". As the angular speed 
co is switched in this embodiment as well, there is no possibility of 
occurrence of a phase-skipping that would be caused if there is a difference 
between the phase calculated at the d-axis detecting unit and the phase 
calculated at the characteristic-feature control unit in switching an angle 9 . 

Next, Fig. 29 is a block diagram, showing the details of the control 
device 1 of Fig. 1, in accordance with the seventh embodiment of the 
invention. 

In common with the first embodiment to the sixth embodiment, the 
voltage vector locus of the high-frequency current to be superimposed draws 
a perfect circle and the d-axis phase angle is corrected on the basis of at 
least one of the length of a long axis and that of a short axis of the current 
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vector locus of the above high-frequency current. While, in the seventh 
embodiment, the current vector locus of the high-frequency current to be 
superimposed draws a perfect circle and the d-axis phase angle is corrected 
on the basis of at least one of the length of a long axis and that of a short 
axis of the voltage vector locus of the above high-frequency current. 

In Fig. 29, a frequency separating unit 24 separates a high-frequency 
current from a three-phase currents inputted from the current sensor 4, by 
means of a general frequency filter. 

A high-frequency current control unit 25 inputs the above 
high-frequency current and a target value in high-frequency current. 
Further, the unit 25 performs PI control etc. so that a high-frequency current 
vector locus draws a perfect circle (see Fig. 30A) to generate a 
high-frequency three-phase voltage command. 

As for the high-frequency control, as mentioned before, it can be carried 
out on the d/q-axes (in the rotating coordinate system) in place of the 
three-phase control (in the static coordinate system). In this way, when the 
high-frequency control is carried out on the d/q-axes while performing the 
detection of a d-axis and the control of a characteristic feature on the 
d/q-axes, there is no slippage of frequency during the rotation of a motor 

When the above high-frequency three-phase voltage command is 
converted to a vector, it follows an elliptic locus as shown in Fig. 30B. In 
a low load zone causing no magnetic saturation, the long axis of the ellipse 
points to a q-axis (i.e. a position of maximum inductance). Thus, the 
detection of a d-axis can be accomplished by deducting 90 degrees from the 
phase of the q-axis. That is, the 3/2-phase converting unit 17 inputs the 
high-frequency three-phase voltage command from the high-frequency 
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current control unit 25 and outputs a voltage in the Of/ /3 -axes of 
coordinates as shown in Fig. 30B. The d-axis detecting unit 28 calculates a 
d-axis phase angle 6 0 by deducting 90 degrees from the phase of the 
q-axis in the above way and outputs the phase angle 9 0 . Noted that the 
above phase angle of the long axis (i.e. q-axis) of the ellipse can be 
calculated by using an expression that results from the replacement of the 
members of I a , I j3 , V a and V /3 in the expression (1) by the members 
of <t> i a , 4> i # , <$> v a and <i> v j3 , respectively. 

As for the correction by a characteristic feature in a high load zone, a 
phase is corrected by maintaining a predetermined value of a characteristic 
feature based on at least one of the length of a long axis of the ellipse of Fig. 
30B and that of a short axis. This characteristic feature is shown in Figs. 9, 
11, 13, 15, 17, 19 and 21. 

Note, in the properties of Figs. 9, 11, 13, 15, 17, 19 and 21, the 
characteristic features have upward-sloping features as the drive current 
phase angle j3 increases. Therefore, if an actual characteristic feature is 
smaller than its target value, the detection phase is advanced (increasing of 
the correction angle 0 ') so as to advance the rotating coordinate system 
while fixing the drive current phase angle j3 at a constant, thereby 
producing an effect to put the current phase angle /3 forwardly. While, if 
an actual characteristic feature is larger than its target value, the detection 
phase is delayed (decreasing of the correction angle 0 ') so as to retard the 
rotating coordinate system while fixing the drive current phase angle j3 at 
a constant, thereby producing an effect to put the current phase angle j3 
backwardly. Then, the control gains are determined by way of experiments. 
Note, as the feature (i.e. the characteristic of the calculation (a/b) of Fig. 15, 
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if a characteristic feature has a feature of downward-sloping as the drive 
current phase angle j3 increases, the above-mentioned relationship 
between advance and delay is reversed. 

Fig. 31 is a flow chart of the process to start the torque control in case of 
no magnetic saturation (out of a high load zone), corresponding to the 
operation of Fig. 4B. Note, the successive steps shown in Fig. 31 follow 
the flow chart of Fig. 4B. 

In the flow chart of Fig. 31, at step S311, a voltage outputted from the 
high-frequency current control unit 25 is converted onto the a / /3 -axes by 
the 3/2-phase converting unit thereby to establish a state as shown in Fig. 3B. 
At sequent step S312, as for the a -axis component v a and the /3 -axis 
component v # of a voltage vector converted in the CL / {$ ~ coordinate 
system, their respective peak values are detected to obtain their amplitudes 
and phases, as shown in Fig. 32. Subsequently, the operations at step S313 
to step S3 17 are similar to those executed at step S7 to step Sll of Fig. 4B. 

In case of detecting the long axis from a voltage vector locus on the a / 
/3 -axes as this embodiment, the expressions (3), (9), (10) and (11) in the 
first and second embodiments will be modified to the following expressions 
(17), (18), (19) and (20), respectively. 

a= {V a 2 +V fi 2 +2V 0 V 5 | sin(T/> Va -^v/3 ) I } 

+ AT{V a 2 +V^ 2 - 2V a V fl |sin(tfva-^vfl ) I }] 

b= {V a 2 +V;+2V a V 6 2 | sin^va-^ve )l} 

-AT{V a 2 +V 5 2 - 2V a V 0 I sin(i/> Va -^v0 ) I }] 
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(17) 



a + b= /" {V a 2 +V fl 2 +2V tt V 5 | sin(^va-^v^ ) I } 

- (18) 

V"(a 2 + b 2 )= V" (V a 2 +V, 2 ) 

••• (19) 

ab = V a V 0 | sin(^ Va -^vfi ) I 

••• (20) 

Next, Fig. 33 is a block diagram, showing the details of the control 
device 1 of Fig. 1, in accordance with the eighth embodiment of the 
invention. 

In the correction of a d-axis phase angle on a basis of at least one of the 
length of a long axis of a voltage vector locus on the a / j3 -axes of a 
high-frequency current and the length of a short axis of the same, according 
to the embodiment, a d-axis phase angle 6 0 is calculated from an a / /3 - 
axes voltage and a phase angle 6 ' is calculated by integrating the detected 
angular speed co ' detected from the characteristic feature". Further, either 
the d-axis phase angle 0 0 or the phase angle 9 ' is outputted in switching. 
The constitution of this embodiment corresponds to that of the second 
embodiment (see Fig. 22) employing a current vector. 

In Fig. 33, the d-axis detecting unit 18 calculates the d-axis phase angle 
0 o based on an a / j3 -axes voltage vector locus. While, at the 
characteristic-feature control unit 20', an angular speed co '[rad/s] is 
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calculated. Then, the so-calculated angular speed co ' is integrated to 
calculate a phase angle 6 ' [rad] by an integrating circuit 21. With the 
exception of high load zone, the switch SW1' is connected with the d-axis 
detecting unit 28 to output the d-axis phase angle 9 0 as a phase angle. In 
the high load zone causing magnetic saturation, the switch SW1 ' is 
connected to the integrating circuit 21 to output the phase angle 6 ' as a 
phase angle. 

Noted that the angular speed oo is calculated by the above expression 
(12). 

Next, Fig. 34 is a block diagram, showing the details of the control 
device 1 of Fig. 1, in accordance with the ninth embodiment of the 
invention. 

In the arrangement to correct a d-axis phase angle based on at least one 
of the length of a long axis of a voltage vector locus on the d/q-axes of a 
high-frequency current and the length of a short axis of the same, according 
to the embodiment, a d-axis phase angle 6 0 is calculated from a d/q-axes 
voltage by employing an angular speed CO . 

The constitution of this embodiment corresponds to that of the third 
embodiment (see Fig. 23) employing a current vector. 

According to the ninth embodiment, the d-axis/q-axis converting unit 27 
converts a high-frequency voltage to a d/q-axes voltage on the d/q-axes 
conversion. While, at the d-axis detecting unit 28', both amplitude and 
phase of the voltage vector are determined by detecting respective peak 
values and zero-cross periods of the d-axis component vd and the q-axis 
component vq of the voltage vector locus converted onto the d/q-axes. At 
this time, although an angle 6 between the d-axis and the long axis of an 
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ellipse is similar to that of the above expression (1) expect to alter the 
members of I a , I (3 , <f> i a and 0 i /3 to the members of vd, vq, (f> vd 
and <f> vq, respectively, it is executed to multiply the angle 9 by a regular 
coefficient (i.e. proportional calculation) thereby outputting an angular 
speed co . Then, in the integrating circuit 22, the angular speed CD is 
integrated to generate a d-axis phase angle 9 0 for output. 

Note, as similar to the third embodiment (Fig. 23), the angular speed co 
may be calculated by a PI (proportional integrating) calculation. Or again, 
on condition of having map-tables with every torque each defining the 
relationship between a torque detected by an experiment in advance and the 
corresponding angle 9 from the d-axis to the long axis of an ellipse and 
further determining the angle 6 corresponding to a target torque by means 
of the specified map-table, it is possible to improve the accuracy of an 
estimated phase. 

The characteristic- feature control unit 20 outputs a correction angle 9 9 
from a d/q-axes voltage. In a high-load zone where the switch SW1 is 
turned ON, the characteristic-feature control unit 20 corrects the d-axis 
phase angle 9 0 by the correction angle 9 ' to provide a phase angle 6 , as 
similar to Fig. 23. 

Fig. 35 is a flow chart of the process to start the torque control in case of 
no magnetic saturation (out of a high load zone), corresponding to the 
operation of Fig. 4B. Note, the successive steps shown in Fig. 35 follow 
the flow chart of Fig. 4A. 

In the flow chart of Fig. 35, at step S351, a voltage outputted from the 
high-frequency current control unit 25 is converted onto the d/q-axes 
thereby to establish a state as shown in Fig. 30B. At sequent step S312, as 
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for the d-axis component vd and the q-axis component vq of a voltage 
vector converted in the d/q-coordinate system, their respective peak values 
are detected to obtain their amplitudes Vd, Vq and phases <t> vd, <f> vq, as 
shown in Fig. 36. Subsequently, the operations at step S353 to step S355 
are similar to those executed at step S243 to step S245 of Fig. 24. 

Similarly to the third embodiment, due to the detection of the long axis 
on the d/q-axes, the present embodiment has an advantage of eliminating 
detection errors since the locus of a high-frequency voltage vector comes to 
a standstill while drawing an ellipse on the d/q-axes when the motor rotates 
at regular revolutions. 

Besides, in detecting the long axis on the d/q-axes as this embodiment, 
the expressions (3), (9), (10) and (11) in the first and second embodiments 
will be modified to the following expressions (21), (22), (23) and (24), 
respectively. 

a = ![/" {V d 2 + V q 2 + 2V d V q | sin (ip Vd - Vq ) I } 

+ ->f{V d 2 +V q 2 - 2V d V q I sin(^vd-^v q ) I }] 
b= {V d 2 + V q 2 +2V d V q 2 I sin(^vd-^v q ) I } 

-•xf{Vd 2 + V q 2 - 2V d V q I sin(^vd-^v q ) I }] 

- (21) 

a + b= -f {V d 2 +V q 2 +2V d V q | sin(^vd-^vq) I } 

- (22) 
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V"(a 2 + b 2 )= /" (V d 2 +V q 2 ) 



(23) 



ab = V d V q | sin (^ yd" ^ v q ) I 

••• (24) 

Next, Fig. 37 is a block diagram, showing the details of the control 
device 1 of Fig. 1, in accordance with the tenth embodiment of the 
invention. 

In the arrangement to correct a d-axis phase angle based on at least one 
of the length of a long axis of a voltage vector locus on the d/q-axes of a 
high-frequency current and the length of a short axis of the same, according 
to the embodiment, a d r axis phase angle 6 0 is calculated from a d/q-axes 
voltage by employing an angular speed co . In addition, a d-axis phase 
angle 6 ' is also calculated from a characteristic feature with an angular 
speed co'. In this embodiment, either the d-axis phase angle 9 0 or the 
d-axis phase angle 6 ' is outputted in switching. 

In Fig. 37, the d-axis detecting unit 28', both amplitude and phase of the 
voltage vector are determined by detecting respective peak values and 
zero-cross periods of the d-axis component vd and the q-axis component vq 
of the voltage vector locus (after separation of high frequency) converted 
onto the d/q-axes. Then, the amplitude and phase are multiplied by a 
regular coefficient thereby outputting an angular speed co . Then, in the 
integrating circuit 22, the angular speed co is integrated to generate a 
d-axis phase angle 9 0 for output. 

On the other hand, the characteristic-feature control unit 20' outputs an 
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angular speed co ' [rad/s] and thereafter, it is further integrated to calculate 
a phase angle 6 ' [rad] by the integrating circuit 21. 

With the exception of high load zone causing magnetic saturation, the 
switch SW1 ' is connected with the integrating circuit 22 to transmit the 
d-axis phase angle 6 0 from the circuit 22 to the succeeding circuit, as the 
phase angle 6 . In the high load zone causing magnetic saturation, the 
switch SW1' is connected to the integrating circuit 21 to transmit the phase 
angle 0 ' obtained as a result of integrating the angular speed co ' to the 
succeeding circuit, as the phase angle 6 . 

Next, Fig. 38 is a block diagram, showing the details of the control 
device 1 of Fig. 1, in accordance with the eleventh embodiment of the 
invention. 

In this embodiment, the integrating circuit 23 for integrating an angular 
speed co in the tenth embodiment is connected to the backside of the 
switch SW1\ 

In Fig. 38, an angular speed co from the d-axis detecting unit 28 5 is 
transmitted to the switch SW1' directly, while the characteristic-feature 
control unit 20' outputs an angular speed co as similar to the embodiment 
of Fig. 37. Then, the switch SW1' is connected with the 
characteristic-feature control unit 20' in a high-load zone causing magnetic 
saturation. Except the high-load zone, the switch SW1' is connected with 
the d-axis detecting unit 28'. Therefore, in the high-load zone, the 
integrating circuit 23 integrates the angular speed co ' from the 
characteristic-feature control unit 20' to output a phase angle 6 . Except 
the high-load zone, the integrating circuit 23 integrates the angular speed co 
from the d-axis detecting unit 28' to output a phase angle 6 . 
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According to the constitution of Fig. 38 5 as similar to the fifth 
embodiment (see Fig. 27), since the switching operation is carried out in the 
form of an angular speed co , it is possible to solve the above problem of 
phase-skipping in switching the phase angle 8 . 

Next, Fig. 39 is a block diagram, showing the details of the control 
device 1 of Fig. 1, in accordance with the twelfth embodiment of the 
invention. 

In the correction of a d-axis phase angle on a basis of at least one of the 
length of a long axis of a voltage vector locus on the a / j3 -axes of a 
high-frequency current and the length of a short axis of the same, according 
to the embodiment, an angular speed co ' is corrected by a correction 
angular speed co The constitution of this embodiment corresponds to a 
constitution resulting from the application of the sixth embodiment (see Fig. 
28) to a "voltage vector" type constitution. Noted that the calculation of 
the correction angular speed co" employs the expressions (12) and (12'), as 
similar to the sixth embodiment. 

In Fig. 39, the angular speed co outputted from the d-axis detecting 
unit 28 ' is transmitted to the accumulator 26 directly. While the 
characteristic-feature control unit 20" outputs the correction angular speed 
co In a high-load zone where the switch SW1' is turned on, the 
correction angular speed co " is added to the angular speed co at the 
accumulator 26. 

The integrating circuit 23 integrates an output from the accumulator 26 
and further outputs a phase angle 0 . Therefore, except the high-load zone, 
the integrating circuit 23 integrates the angular speed co outputted from the 
d-axis detecting unit 28 ' to output the phase angle 9 . While, in the 
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high-load zone, the integrating circuit 23 integrates a value obtained by 
correcting the angular speed co by the correction angular speed co " to 
output a phase angle 6 . Also in this embodiment, as the angular speed 
co is switched, there is no possibility of occurrence of a phase-skipping that 
would be caused if there is a difference between the phase calculated at the 
d-axis detecting unit and the phase calculated at the characteristic-feature 
control unit in switching an angle 9 . 

Finally, it will be understood by those skilled in the art that the foregoing 
descriptions are nothing but some embodiments of the disclosed control 
device for an electric motor. Besides these embodiments, various changes 
and modifications may be made to the present invention without departing 
from the spirit and scope of the invention. 

Japanese Patent Applications No. 2002-222623 filed on July 31, 2002 
and No. 2003-271330 filed on July 7, 2003 are incorporated by reference 
herein in its entirety. 

While the present invention has been disclosed in terms of the preferred 
embodiments in order to facilitate better understanding thereof, it should be 
appreciated that the invention can be embodied in various ways without 
departing from the principle of the invention. Therefore, the invention 
should be understood to include all possible embodiments and modifications 
to the shown embodiments which can be embodied without departing from 
the invention as set forth in the appended claims. 

The scope of the invention is defined with reference to the following 
claims. 
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